Electroweak jet cascading in the decay of superheavy particles 
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We study decays of superheavy particles X into leptons. We show that they initiate 
cascades similar to QCD parton jets, if mx > 10 6 GeV. Electroweak (EW) cascading is 
studied and the energy spectra of the produced leptons are calculated in the framework 
of a broken SU(2) model of weak interactions. As application, important for the Z- 
burst model for ultrahigh energy cosmic rays, we consider decays of superheavy particles 
coupled on tree-level only to neutrinos and derive stringent limit for these decays from 
the observed diffuse extragalactic 7-ray flux. 
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A characteristic feature of high energy deep- 
inelastic scattering, e + e~-annihilation and decays 
of superheavy particles is the cascading of QCD 
partons. Despite the smallness of the QCD cou- 
pling as, cascading occurs because the probabil- 
ity of parton splitting is enhanced by large loga- 
rithms for soft or collinear parton emission. Re- 
cently, similar logarithms were found to dominate 
also the EW radiative corrections at the TeV scale 
and above 0. Evolution equations for these cor- 
rections, similar to the DGLAP equations in 
QCD but valid for a spontaneously broken theory, 
were derived in Q. In this Letter, we want to 
draw attention to decays of particles much heav- 
ier than the EW scale: if these particles decay 
or annihilate, then the effects discussed above re- 
sult in particle cascades which proceed through 
all U(l) x SU(2) x SU(3) interactions with couplings 
evi, ct2 and 03 enhanced by logarithmic terms. 
Particular attention is given to EW Lepton-Boson 
(EWLB) cascades developing without QCD par- 
tons. 

We shall illustrate the development of EW cas- 
cades by the example of a superheavy V-particlc 
with mass mx < TOgut decaying into leptons, 
which has an interesting physical application to 
the Z-burst model (see below). Let us discuss 
first the case when the X-particle has on tree- 
level only couplings to lepton pairs 11 and con- 
sider the decay mode X — > vv. For mx 3> rnz, 
the mass of the Z boson is negligible compared to 
the available momentum flow Q 2 < m 2 x /A. Then, 
soft and collinear singularities generate large log- 
arithms In (m^/mfj), which can compensate the 
smallness of a 2 /(27r). 



To show the existence of the cascade, we con- 
sider the ratio R = T(X -> vvZ)/T(X -> vv). 
Neglecting terms finite in the limit mz — > 0, it is 
given by 

R = Z^T- (ln 2 e + 31ne + ...) , (1) 
° 7TC w 

where e = (mz /mx) 2 and c 2 ^ = cos 2, dw- 
The In 2 e term — which compensates the small 
coupling — is typical for soft and collinear singu- 
larities in the emission of vector particles g. For 
mx ~ 10 6 GeV, the decay probability into vv and 
Zvv becomes comparable, R ~ 0.5, signaling the 
break-down of perturbation theory: the decay of a 
particle with mass mx 3> 10 6 GeV, even if coupled 
only to leptons, initiates a cascade, very similar to 
that known in QCD. 

For the quantitative study of EWLB cascades, 
we consider a broken SU(2) gauge theory as a sim- 
plified model for the electroweak sector of the Stan- 
dard Model (SM). Its particle content is a triplet 
of gauge bosons Wj? with mass mw, a physical 
Higgs h with mass rrih and n g — 3 generations of 
leptons / = (Il,W,l)> Ir, with I = e, /j,, r. The 
SU(2) doublet / is coupled to the gauge sector as 

92tthYW-i")fcW*. 

Since the terms generating the logarithms in 
Eq. (Q) are associated with collinear and/or soft 
emission of additional particles, we introduce split- 
ting functions Pijk(z), viz. 

dT n+ i a 2 dt 
dr ~ aui^jk = ——dzP ljk (z) , (2) 

where dT n is the decay width into n particles 
with virtuality t, z is the energy fraction of 
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the additional particle, and dcoi^jk approximates 
dT n+ i/dT n in the high-energy, small-angle scatter- 
ing limit. 

In the following, it is useful to distinguish be- 
tween transverse and longitudinal modes of the 
gauge bosons, Wt and Wl- Going from QCD 
to our broken SU(2) model, there are only trivial 
changes for the transverse modes due to the differ- 
ent structure constants and the V-A coupling, 



P 



WtWtWt 



= 2 



1 - z 



+ z(l-z) 



p WTf f='-f[z 2 +(i-zn 



ffw T 



2 P ff'w± 



3P 



3 l + z 2 



ffw° 



4 1 - z 



(3) 
(4) 
(5) 



External longitudinal modes Wl can be replaced 
by the corresponding Goldstone bosons xw in S- 
matrix elements for s S> myy Since mw acts 
as an infrared cut-off and is assumed to be well 
above the experimental resolution, virtual correc- 
tions and real bremsstrahlung do not mix. Re- 
placing Wl by the scalars xw makes clear that 
the splitting functions involving Wl do not con- 
tain terms singular in z. The longitudinal modes 
Wl are therefore subdominant compared to Wt 
and we will neglect them. 

Finally, we discuss whether the Higgs partici- 
pates actively in the cascade evolution. The Higgs 
can split into pairs of fermions or gauge bosons. 
In the first case, all Yukawa couplings except the 
one to the top quark t can be neglected; moreover 
collinear singularities are absent. Therefore, only 
the processes t — > th and h — > tt have a minor im- 
portance at Q 2 close to the GUT scale. In the sec- 
ond case, i.e. when splittings involve Higgses and 
gauge bosons, we cannot define splitting functions 
in the conventional sense. Instead, dojh^w L w L and 
duh^w T w T are proportional to g 2 /t 2 . The result- 
ing branching probability is negligible. 

Hence, the EW cascade is supported essentially 
by the same splitting modes as in the unbroken the- 
ory. The proof of the factorization of both collinear 
and soft singularities in the EW theory follows the 
same logic as in QCD [|]j3| and the effects of the 
spontaneous symmetry breaking result at s 3> myy 
only in subleading corrections Q. Thus one ob- 
tains the EW evolution equations in the Altarelli- 
Parisi form [Q and is able to establish the SU(2) 
charge coherence picture in the usual way Q . The 
evolution equations in integral form including co- 
herence effects are 

(Q 2 , x, t) = A l s (Q 2 , t) 51 5(1 -x)+ (6) 
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xa 2 [z 2 (l - z) 2 t'] P iM (z) D k ^(x/z,z 2 t') , 

where Di^j(Q 2 ,x,t) defines the energy distribu- 
tion of particles j with energy fraction x at scale 
t produced by the parent particle i at scale Q 2 . 
Here, £i(t') ~ mi/y/i, where mi is the mass of 
the particle I, t = q 2 /[z(l — z)], q 2 the vir- 
tually of the parent particle and z the energy 
fraction of the produced one. We use aa(t) = 
a 2 (m%)/[l + (6/47r)ln.(t/m|)a 2 (m|)] with b = 
19/6 for one Higgs doublet and 0.2(111%) ~ 1/29.6. 

In Eq. (^J), we have introduced the Sudakov 
form factors A\(Q 2 , t) giving the probability of no 
branching for the parton of type i in the scale range 
between Q 2 and t, 
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x ^a 2 [z 2 (l-z) 2 t']P lkl (z). 

Equations (^0) are very similar to the analo- 
gous ones in QCD and allow to use the probabilis- 
tic scheme of Ref . fl . There is however an impor- 
tant difference between the EWLB and QCD cas- 
cades in the physical meaning of the cut-off value 
Ek(t')- In QCD, it is defined by the scale g^ in at 
which the perturbative evolution of the cascade is 
terminated, £qcd(0 = \/?minA'- % contrast, 
q 2 lin of the EWLB cascade is given by the physi- 
cal masses of the particles produced in the split- 
ting i -> kl. For W T -> //, e k (t') w 0, i.e., it 
is much smaller than the one in Wt — * WtWt, 
where e k (t') ~ m\y/VV- As a consequence, the 
Wt does not dominate the cascade evolution as 
strongly as gluons in QCD and the number of Wt 
leaving the cascade is smaller than the number of 
lcptons. 

The EWLB cascade does not terminate abruptly 
at a certain Q 2 . Instead, the probability for 
no further branching increases smoothly from 
Ai(Q 2 ,0) = for Q 2 -► 00 to Ai(0,0) = 1. Never- 
theless, one can define Q 2 T by A l s (Q 2 , 0) = 0.5, i.e. 
by that virtuality Q 2 r at which the probabilities 
for further branching and for non-branching be- 
come equal. Then Q CI is given by Q cr ~ 10 6 GeV 
and the probability for further splittings decreases 
indeed drastically for lower Q. 

We choose to solve the evolution equation with a 
Monte Carlo simulation. This method gives us the 
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advantage of including non-abelian charge coher- 
ence effects through angular ordering. In Fig. [j], 
we show the spectra dN/dl, where I = ln(l/x) and 
x = 2E/mx, of leptons and W's for rax = 10 10 , 
10 13 and 10 16 GeV and for the decay mode X -> 
vv —* all. The average number of produced lep- 
tons and W's increases from 9 for rax — 10 10 GeV 
to 47 for rax = 10 16 GeV, where we consider the 
/i ± , t~^~ and the W a as stable particles. The first 
bin contains between 7% (rax = 10 10 GeV) and 
2%o (mx = 10 16 GeV) of prompt neutrinos with 
E = mx /2 from no-cascading decays. Leptons 
which stop branching after the first splitting pro- 
duce a tail superimposed to the usual Gaussian, 
clearly visible for I < 3. 




J = ln(l/x) 



FIG. 1. Energy spectra dN/dl of leptons / 
(solid lines) and W's (broken lines) produced in 
X -> vv -> all; for m x = 10 16 , 10 13 and 10 10 GeV 
from top to down. 

Let us discuss now the inclusion of QCD partons 
in our model. Since the EW gauge bosons split 
also into quarks q, there will be a mutual trans- 
mutation of "leptonic" and "QCD" cascades. The 
shape of the hadron energy spectra is however only 
marginally influenced by the leptons, first because 
the QCD cascade is determined mainly by gluons 
g and second because the probability of q — ► q + g 
is much larger than of q — > q + W. On the other 
hand, the splittings Wt — > qq act continuously as 
a sink for the particles and energy of the EWLB 
cascade. These splittings are taken into account in 
the MC simulation and the spectra shown in Fig. 1. 

We have calculated also (see Table 1) a quantity 
interesting for decays of V-particles in the Uni- 
verse: the fraction of energy / cm , transferred to 
electrons and photons in the decay mode X — > vv 



as function of mx ■ These particles initiate electro- 
magnetic (e-m) cascades in the intergalactic space, 
through interactions with microwave (CMB) pho- 
tons. In the calculation of f em , we assumed that 
all hadrons produced in Wt — > qq — > hadrons are 
pions, while we used the SM branching ratios for 
/i and r decays. 

The above calculations have interesting con- 
sequences for Ultra High Energy Cosmic Rays 
(UHECR) produced by superheavy Dark Matter 
(DM) (H} and, in particular, by Z-bursts 0. In 
the latter model, UHECR are produced through 
the resonant production of Z-bosons in the colli- 
sions of UHE neutrinos with DM neutrinos, v + 
v — > Z° . The resonant energy of UHE neutrino is 
E = m 2 z /2m v — 4.2 x 10 21 m~^ eV, where m c v is 
the mass of the DM neutrino in eV. 

Decays of these Z-bosons produce UHECR 
which energy spectrum, according to recent cal- 
culations |lj,[l3|], has a weak Greisen-Zatsepin- 
Kuzmin cutoff and explains well the observations. 
A remarkable feature of this model is that the pro- 
duction of UHE particles with energies higher than 
10 20 eV does not involve exotic elementary particle 
physics, while its drawback consists in the neces- 
sity of an enormous flux of resonant neutrinos. 

A widely discussed possibility is the production 
of UHE neutrinos in astrophysical sources. This is 
an unrealistic option because it requires too high 
luminosities of the astrophysical sources: we can 
estimate the required flux of resonant neutrinos 
I v (Eq) and therefrom the neutrino energy density 
u v as uj v « (2.4 - 3.6) x KT 13 /?^ 5 erg/cm 3 . In 
these calculations, we have used that for the flat 
spectra generated by heavy particle decays, such 
as Z-bosons and V-particles, UHE photons domi- 
nate at the highest energies E > 1 X 10 20 eV 
The resulting neutrino luminosity of a source, es- 
timated as L v ~ u) v / (n s to), where n s is the source 
density and to the age of the Universe, is too high: 
(8 — 12) x 10 44 erg/s, if the sources are normal 
galaxies, and (8 — 12) x 10 46 erg/s in the case of 
Seyfert galaxies. Our result confirms similar con- 
clusions of Ref . (r| , reached on the basis of a limit 
due to the diffuse 7-ray background. 

Alternatively, Z-burst neutrinos can be gener- 
ated by decays of superheavy particles, either ex- 
isting as DM particles or generated by topological 
defects. In the usual models this possibility is ruled 
out by the diffuse 7-ray background |l6| . However, 
in Ref. [O] it is suggested that this limit can be 
evaded, if X-particles decay exclusively to neutri- 
nos, i.e. X — > vv. 
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We shall demonstrate now that electroweak cas- 
cading rules out this last hope, too. 

EW cascading modifies the e-m cascading up- 
per limit for UHE neutrinos |l6],[L7| by the fac- 
tor /j/// em , where f v is the fraction of energy mx 
transferred to neutrinos. Then the upper limit on 
the UHE neutrino flux is given as 



47T/ cm .E 2 



(8) 



with w cas w 2 x 1CT 6 eV/cm 3 according to EGRET 
observations. Taking I v (Eq) from Eq. (||), one 
can calculate the flux of UHE photons produced 
in Z-bursts. At E > 1 x 10 20 eV, this flux 
should be of the order of the observed UHECR 
flux. With / cm = 0.2 from Table 1, we obtain 
E 3 I 1 (E) « (5 - 8) x 10 21 m c v eV 2 /m 2 s sr, which 
is almost three orders of magnitude less than ob- 
served. 

Summary — EW cascading in decays of super- 
heavy particles results in EWLB cascades similar 
to the cascade of QCD partons, if mx > 10 6 GeV. 
Thereby the flux of prompt neutrinos, e.g., in an- 
nihilations of DM particles is reduced. EW cas- 
cades allow a probabilistic interpretation and ex- 
hibit destructive coherence at small x. The gen- 
eration of electrons and photons, which are able 
to start e-m cascades on the CMB, in the EW 
cascade allows to exclude those Z-burst models, 
in which the resonant neutrinos are produced in 
X — ► vv decays. Thus, the production of resonant 
neutrinos in astrophysical sources and by decays 
of superheavy particles (both as DM particles or 
produced by topological defects) result in neutrino 
fluxes too low for the Z-burst model. The only 
possibility left is the oscillation of sterile neutri- 
nos, e.g. in hidden-sector/mirror models Jl6| |, into 
ordinary ones. 
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TABLE I. Energy fraction / cm transferred to elec- 
tron and photons in the decay X — > vv of a particle 
with mass mx- 
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